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Abstract 

An improved weighting algorithm apphed to hadron showers has been developed 
for a fine grained LAr calorimeter. The new method uses tabulated weights 
which depend on the density of energy deposited in individual cells and in a 
surrounding cone whose symmetry axis connects the interaction vertex with the 
highest energy cluster in the shower induced by a hadron. The weighting of 
the visible energy and the correction for losses due to noise cuts are applied in 
separate steps. In contrast to standard weighting procedures the new algorithm 
allows to reconstruct the total energy as well as the spatial energy deposition on 
the level of individual calorimeter cells. 

The linearity and the energy resolution of the pion signal in the momentum inter- 
val 2 GeV/ c < p < 20 GeV/ c studied in this analysis are considerably improved in 
comparison to the standard weighting algorithm as practiced presently by the HI 
collaboration. Moreover the energy spectra reconstructed with the new method 
follow in a broad interval a Gaussian distribution and have less pronounced tails. 
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1 Introduction 



The different response of hadron calorimeters to electromagnetic and hadronic 
showers, observed for the first generation of hadron calorimeters 12 El, are a 
source of non-linearities, deviations of the relative resolution from the statistically 
expected I/a/E -behaviour and the origin of the non-gaussian shape of the re- 
sponse function |4, 5j. Two different methods have been proposed to avoid these 
disadvantages. Hardware compensation, optimized by extensive simulations of 
hadron showers [SlElEl, enhances the hadron signal by detecting neutrons, pro- 
duced in hadron ^^^JJ interactions, with scintillators jHlini- On the other hand 
in fine grained calorimeters the signal of different detector cells can be weighted 
in such a way that the overall signal of electrons and hadrons depositing the 
same energy is equalized. This method was pioneered by the CDHS Collab- 
oration |3] using a large scintillator-iron calorimeter to detect hadron showers 
produced in t/-events. Essential improvements were achieved by the HI Collab- 
oration pm [TT| IT^ IT^ exploiting the high granularity of its 45000 channel LAr 
calorimeter. 

In this paper a new weighting algorithm is presented which not only allows to 
determine the deposited energy in an improved way but in addition permits a 
more realistic reconstruction of the spatial shower distribution. While previous 
algorithms aimed just for a determination of the total shower energy, the new 
method intends to reconstruct the energy in each cell of the calorimeter hit by the 
hadronic shower. Losses due to noise suppression cuts applied to individual cells 
are considered separately by adding a correction term to the total reconstructed 
energy. 

After a short description how the visible energy is derived from the measured 
charge, which is identical to the standard procedure of the HI algorithm ^| ^1 
IT7j . the new method is explained. The weights, and the measurable variables 
they depend on, are determined next, followed by a description of the algorithm, 
correcting for losses due to noise. Finally the new weighting algorithm is applied 
to Monte Carlo and test beam data. The results achieved are compared to those 
of the standard HI algorithm described in [TH UHl HZ] . 
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2 Reconstruction of the Visible Energy 



2.1 Determination of the Electromagnetic Calibration 
Constant 

The electromagnetic calibration constant which allows to convert the recorded 
charge into visible energy, has been derived by an iterative procedure ^] equal- 
izing the reconstructed energy of test beam and simulated data for electrons: 

{E^) = {Er) (1) 

The experimentally determined energy is derived from the deposited charges Qi 

(Br) = J2 (2) 

i 

where the sum runs over all calorimeter cells with a signal passing the noise cuts. 
The reconstructed signal of simulated data is given by the expression 

{Eq"^) = ( ^ {Cs^m El^^ + Cexp Qnoise) ) (3) 
i 

C„™ = {^) (4) 

Csim is derived from simulated data. Edep is the energy deposited in the calorime- 
ter, while 

= J2 K^s (5) 

i 

is the detected visible energy in the active medium of the calorimeter derived by 
simulation. The influence of saturation phenomena, charge losses due to electron 
capture by electronegative gases |^, recombination in the ionization column [20] 
etc. was derived from HV curves recorded before and after the data collection. In 
addition probes similar to those used finally in the HI calorimeter monitored 
the response to the signal of a radioactive source continuously. Qnoise is the 
measured noise in calorimeter channel i. Cexp is obtained iteratively from eq. (j21) 
- requiring the constraint to hold for electrons. 
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The calibration constant Cexp determined by this method defines the electro- 
magnetic scale of the calorimeter. The reconstructed energy of cell i on the 
electromagnetic scale for test beam data is given by the expression 



Eq = Cexp Qi (6) 

and for simulated data by 



Eq — Csim Eyj^g + Cexp Q noise i'^) 

in analogy to eqs. (jSj) and 0- 

The data used in the present analysis were collected at the CERN SPS test 
beam H6 |22] using a FB2-type module of the HI LAr calorimeter 121] • The 
simulation employed the HI software packages HISIM [23] and ARCET [211 based 
on the event generators GEANT 3.21 [2^ and GHEISHA[2Hj. 



2.2 Clustering of Energy 

The identification of calorimeter cells hit by the hadron showers follows the stan- 
dard methods developed by the HI Collaboration, described in more detail in ref. 
[m CHI CZI • Besides the charge produced by the hadron shower electronic noise 
contributes to the detected signal which amounts typically to anoise = 15 MeV 
up to 30 MeV per channel Only charges with 

U\ > 2.5 (8) 

are recorded for the calorimeter modules used in the present analysis. In the 
analysis however an additional cut is applied. Besides those channels with 



m > ^<oise (9) 

in addition cells are taken into account which are direct neighbours of a channel 
withQi > 4cT;„i,^: 



Qi> A (jj A \Qj\ > 2.5 (Tj with j = neighbour cell of i (10) 

Clusters combining those cells which pass the noise cuts (jS)), ffTUI) are con- 
structed in two steps ^3]. First all cells in a given plane of the calorimeter 
at constant distance from the beam are grouped into two-dimensional clusters 
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around the cell with the highest charge deposit. In the second step adjacent two- 
dimensional clusters are combined to a three-dimensional cluster. The topological 
nearest cells to a 3D-cluster with a negative signal are added to the original clus- 
ter. If the total charge of this combination is negative the 3D-cluster is excluded 
from the following analysis. 

Clusters due to primary photons and electrons are identified by estimators which 
exploit the characteristic shape of an electromagnetic shower fT^ I16j. They are 
not considered further in the analysis. The remaining clusters are labeled as 
hadronic clusters, all weights are derived from and applied to this subset of the 
data. 

The quantity Econe, the energy deposited inside a cone whose axis connects the 
interaction vertex with the most significant hadronic cluster, where its significance 
is defined by 




turns out to be an important input variable for the weighting algorithm developed 
in this paper. Econe is determined in two steps following essentially the procedure 
developed for the HI standard analysis [in]. If no other cluster apart from the 
one defining the cone axis is found within an opening angle of 11° the energy of 
the cone is given by 



Econe ^cluster (l^) 

However, if more than one cluster is recorded within this cone, a new one of 11° 
opening angle is constructed whose axis connects the interaction vertex with the 
energy center of gravity of the hadron clusters of the previous cone. In this case 
the cone energy is defined by the sum over all clusters in the new cone 



Econe — ^duster (13) 

j 

In the following the cone reconstructed with the largest energy is referred to as 
"most energetic cone"; E^'^^ denotes its energy. Further cones are constructed 
from the remaining hadron clusters following the procedure described above. 
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2.3 Description of the New Weighting Algorithm 

Two components of the hadron showers strongly influence the energy deposition 
and the possibilities to measure its energy. The break up of nuclei strongly reduces 
the detectable energy [310111; the lost energy is referred to as "invisible energy" 
in the following, it has to be compensated by weighting. On the other hand the 
electromagnetic component in a hadron shower is deposited in the calorimeter 
without losses and therefore has not to be weighted. 

The energy density 



in cell i with volume Vol'^ allows to tag these components jl]. 

In figs. b the fractional contribution of the invisible energy and the electro- 
magnetic energy to the total deposited energy in cell i respectively are plotted 
as functions of the energy density pi. The electromagnetic component rises with 
increasing energy density. This behaviour can be traced back to the fact that 
electromagnetic subshowers have a smaller spatial extension than hadron induced 
subshowers. 

In figO^ the theoretical weights 



<,o = (15) 

which according to eq.Q do not consider the noise contribution are shown as a 
function of pj. It reveals a strong variation with pi resulting from the individual 
dependencies of the invisible and the electromagnetic component on the energy 
density, wlf^^^ is compared in figlHi to the weight 



< = % (16) 

which allows to convert the reconstructed energy of eq.Q on the electromag- 
netic scale into the real deposited energy -E^gp. The differences and wj;^ „ at 
small energy densites pi can be explained qualitatively as a consequence of the 
overlayed noise and the noise cuts applied pUj . 

The theoretical weights wlf^ ^ in the interval 3 < Pi < 30 decrease (figHt;), 
since in this region the relative contribution of nuclear binding energy diminishes 
(fig^), i.e. less energy has to be corrected for. In the interval 0.2 < Pi < 
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3 the theoretical weights wj^^ increase with pi, this can be attributed to 
the growth of the invisible energy (fig^). Especially at lower values of pi the 
energy of charged particles is deposited by excitation and ionization of atoms; 
hence contributions from nuclear reactions can be neglected. In the interval 
Pi < 0.2 ^ the weight wlf^g increases with decreasing pj. Since these energy 
densities are characteristic for energy depositions of low energy electrons and 
photons in the tail of the shower (see increase of the fractional contribution of 
the electromagnetic energy in figUlD), the transition effect [3 El EI] reduces the 
signal and therefore forces the weights to increase. 

Fig|21 demonstrates that in first approximation the shape of the energy density 
distribution is independent of the energy of the primary pion, but a closer look 
reveals differences especially at large values of pi. This can be attributed to 
the increase of the electromagnetic fraction in hadron showers at higher primary 
energies [21]. This effect has to be taken into account in a global way; the energy 
Econe of the cone provides the necessary information to estimate the initial hadron 
energy jSUj as will be demonstrated. 

Hence the weighting factors of the new algorithm are parameterized as function 

of Pi and Econe 



{ED 



(17) 



pi ■jEcon 



They are derived from simulated data taking noise into account and applying 
noise cuts. (E^i^p) and (Eq) are the mean values of the deposited and recon- 
structed energy in a pj — Econe interval respectively. Only cells of the "most 
energetic cone" of an event are considered, since E^'^^ is a convenient measure of 
the hadron energy initiating the shower. This becomes evident from fig|21 where 
the fraction of energy deposited in the "most energetic cone" is plotted as a func- 
tion of the primary hadron energy. For > 2 GeV/c more than 75% of the total 
reconstructable energy is detected in this cone, this fraction rises strongly with 
increasing and saturates at ~ 95%. Hence for the construction of the weights 
the choice of E^^^ as the second variable, on which if^g^ depends, makes sense. 
If the weights of cells for cases with Econe < E^^c would have been considered 
in addition, corrections for showers with low and high deposited energy would 
be mixed up, while considering only the "most energetic cone" provides clean 
conditions with small overlap of low and high energy data (figlH). 

The weights according to equation (|17|) of the new algorithm for the electro- 
magnetic and hadronic modules of the Hl-LAr calorimeter in the test set-up are 
collected in figs. 5a, b. A double logarithmic scale has been chosen to depict the 
strong variation of w^etu small pi, Econe EOl- The new weighting algorithm 
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Figure 1: Fraction of energy deposited in cell i normalized to the corresponding 
reconstructed energy as function of the energy density pf a) invisible energy, 
b) electromagnetic energy. Theoretical weights wlf^^ ^ according to eq. (fT^ are 
plotted as a function of pi in c) and are compared to the weights wlf^ (full points) 
of eq. (fT^ in d). The energy of the primary hadron is 15 GeV. 

allocates for each cell i according to its energy density and the energy of its as- 
signed cone the weighting factor w^g^ (pi, Econe) from the tables shown in figs. 5. 
It reconstructs the energy on the hadronic scale E^.^^ by multiplying the energy 
on the electromagnetic scale E^ with the new weight 



■^rec -^0 ' '^new {Pij Econe) 



Note that the reconstruction of the energy on the hadronic scale of cells, which 
do not belong to the most energetic cone, is also performed with the weighting 
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Figure 2: Distribution of energy density of cells in the "most energetic cone" for 
pions of a) 6 GeV/c and b) 25 GeV/c. 
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Figure 3: Fraction of energy reconstructed in the "most energetic cone" normal- 
ized to the total reconstructed energy. 
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Figure 4: Distribution of reconstructed energy in the "most energetic cone" a) 
for low and b) for high momentum pions. 

factor shown in fig. 5. 

The apphcation of the new weighting algorithm allows to reconstruct the de- 
posited energy on the cell level. This is demonstrated by figEl where the the- 
oretical weights wlf^ of feqlT^ are compared to the applied weights of the HI 
standard and of the new algorithm as described in this paper. While the theo- 
retical (full circles) and the applied weights of the new algorithm (open triangle) 
coincide in good approximation for all values of pi, the HI standard weights 
(open circles) differ in their trend from the expectation. Note however that by 
construction the total energy deposited in the calorimeter is successfully recon- 
structed by both algorithms. The standard HI weighting procedure achieves this 
through the iterative application of weighting functions according to equation: 
-E*gc = Ciexp{—C2Y^) + C'Aemc/hac ' where Cj are parameter functions, 
which are depending on the jet energy and polar angle of the jet |T2]. These 
weighting and parameter functions have been determined in an iterative process 
optimizing both the reconstruction of the total energy deposited in the hadronic 
shower and the resolution. In this iterative process no emphazis is put on the 
reconstruction of the deposited energy on the level of single cells. 
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Figure 5: Weighting factors in percent a) for the electromagnetic and b) the 
hadronic FB2-modules of the HI LAr test calorimeter set-up j22 cis function of 
he energy in the "most energetic cone" and the energy density in cell i. 
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Figure 6: Comparison of the theoretical weights (full circle) according to eqlTHl 
with the new (open triangle) and the HI standard weights (open circle) as func- 
tion of the energy density in a cell for 15 GeV pions. 
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2.4 Correction of Energy Losses due to Noise Cuts 



Only those cells can be considered by the weighting algorithm which have a finite 
signal after the noise cuts, i.e the new weights derived only optimize the energy 
reconstruction for cells with a signal above the noise cuts. Hence the losses due 
to this cut have to be taken into account in a separate step. In contrast to this 
procedure in the HI standard analysis the noise correction is included in the 
weighting factor. 

In the new algorithm a first order approximation noise correction is derived using 
information of cells belonging to the "most energetic cone". All cells of this 
cone are considered including those which do not belong to a hadronic cluster. 
From simulated data including detector noise the total energy deposited in the 
cone, E^^P'^"^, and the deposited energy of cells passing the noise cut E^^^'^^"^ are 
calculated. The correction is given by the expression 

A rpdep,cor rpdep,tot rpdep,rem ^"1 

cone cone cone ' \^) 

The energy correction is evaluated function of 

Efone = E ^rec (20) 

i 

on the electromagnetic scale, where the sum runs over all cells in the "most 
energetic cone". It is added to the weighted energy of the cone. The result is 
shown in figE The weak energy dependence for Ef^^^ < 1 GeV is due to the fact 
that in cones with small energy the fractional noise contribution is large hence the 
signal loss due to the noise cuts is small. With increasing -EcoL; i-^- increasing 
shower energy, more cells in the cone have a signal which can be suppressed by 
the noise cut, hence the infiuence of the noise cuts grows. 



3 Comparison with Simulation and Test Beam 
Data 

The new algorithm is applied to a set of simulated data. The geometry of the 
simulation corresponds to the CERN test beam configuration [22] with an impact 
angle 6 = 33.73° of the primary pion. Also the geometry of the detector and the 
cryostat are adapted to the CERN test set-up [5^. 

Typical distributions of the reconstructed energy are shown in fig|H^ for the stan- 
dard HI and the new improved algorithm (figlHb)- For both weighting procedures 
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Figure 7: Correction defined in eq.(fT^ considering the losses due to the noise 
cuts as a function of the sum of all cells in the "most energetic cone". 

a Gaussian shape of the energy distribution is achieved in the peak region. For 
the new algorithm the tails are smaller, especially at higher beam energies the 
high energy tail is strongly suppressed. Note that energy losses due to material 
in front of the calorimeter are not corrected for, hence the reconstructed and the 
primary energy differ. 

The energy response is linear within 2% for energies above 2 GeV if one uses the 
new algorithm while for the standard HI procedure larger deviation from linearity 
are observed (figE))- An even better linearity, especially in the region of 5GeV, 
would be achievable when using a finer grid in the Monte Carlo simulations for 
the determination of the new weighting factors. The use of the new algorithm 
furthermore leads to an improved energy resolution of the calorimeter ffiglTUI). 

Finally the new algorithm also was applied to real data taken at the CERN SPS 
test beam H6 with a pion beam energy of 20 GeV. In fig^2 the reconstructed 
energy distribution for the two algorithms are compared, again showing an im- 
provement when using the new weigths confirming the results from the study of 
simulated data. 



4 Conclusion 

A new weighting algorithm has been developed to correct the different response of 
the HI LAr calorimeter to electromagnetic and hadronic showers. The essential 
difference between the new and the standard HI algorithm consists in the separate 
treatment of weight factors and the corrections for noise cuts. This separation 
allows to reconstruct the energy deposited in the cells of the calorimeter properly. 
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Figure 8: Comparison of the reconstructed energy distributions a) applying the 
standard HI (left row) and b) the new weights (right row) for vr-mesons of two 
energies of the primary pion (Monte Carlo simulation). 

Hence the new algorithm described in this paper not only improves the linearity 
and the energy resolution (figsM fTUI) but in contrast to the HI standard algorithm 
it allows to reproduce the shower shape (figEl)- 

Recently this algorithm has been generalized by J.Marks |22] in such a way, that 
it can be applied to all modules of the HI LAr calorimeter at the HERA ep 
storage ring. Moreover improved corrections for the noise cuts were developed. 



14 



-g 0.4 
LU 



Q. 
Q) 
T3 

LU 

I 

o 

Q) 
LU" 



0.2 r 







J 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 


1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 




1 1 1 1 1 1 1 




























A 


A 






JILL 





o 


O 


ti — 




JUL 


r 












-1 1 H 1 1 1 1 1 1 


1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 F 





5 


10 


15 


20 





Pj^ [GeV/c] 



Figure 9: Difference between reconstructed and deposited energy normalized to 
the deposited energy as function of tfie momentum of the primary pion for the HI 
standard (open circle) and the new weights (triangle) (Monte Carlo simulation). 
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Figure 10: Relative energy resolution of the reconstructed energies for HI stan- 
dard (circles) and for new weights (triangles) (Monte Carlo simulation). 
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Figure 11: Distribution of reconstructed energies for test beam data collected at 
the CERN test beam H6, data were collected with test modules of the HI LAr 
calorimeter: reconstructed a) with new and b) with HI standard weights. 
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